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ASSTRACT 

Pile oscillators have heen successfully employed to detect 
small concentrations of inpurities in substances such 78s boron in 
graphite. In several instances the technique has proved more 
sensitive than the commonly used "dancer coefficient" method. The 
danger coefficient method using the Aerojet General Nucleonics 201 
Reactor has been demonstrated to be as sensitive as the pile 
oscillator wethod on other reactors; thus a pile oscillator for the 
AGN 201 Reactor offered the possibility of a further improvement. 

Because of low cost and simplicity, a design utilizing one of 
the installed ion chambers in the reactor and the linear amplifier 
on the control console was chosen for construction of the pile 
oscillator. Additional instrumentation employing the output of the 
linear amplifier gave excellent performance in reproducing a test 
siynal; however, when the pile oscillator was put into operation 
satisfactory reproduction of the output for a given sample could not 
be obtained. It is felt that the error is due to reactor noise. 

It is concluded Chater pile oscillator utilizing the installed 
ion cha:..ber and linear amplifier of the AGN 2C1] Reactor operating at 
its present power level is not satisfactory because of reactor noise. 
It is doubtful that any substantial improvement could be obtained 
unless the operating power were increased above the vresently 


licensed value. 
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1. Introduction. 

The time behavior of a nuclear reactor way be affected to 
varying degrees by tue introduction of a material into the reactor, 
This behavior can be related to the nuclear cross eecticee: of the 
material. Methods have been devised to utilize the time behavior 
characteristics of tyne reactor to measure, in particular, neutron 
absorption cross sections. 

Certain important impurities have microscopic absorption cross 
sections which are very large. S:iall percentages in materials will 
increase the total macroscopic cross section significantly. Boron is 
a common exm.ple and its presence in graphite or iron must often be 
quantitatively determined. The methods by which absorption cross 
section measurements are made may readily be adapted for measuring 
the concentration of impurities. This is usually done »y a comparison 
with a standard... 

One technicue is called the danger coefficient method {1} ©. A 
sample of material is placed into the reactor in a fixed position, 
usually the center of the core. The time behavior of the reactor is 
noted and through appropriate relationships a comparison made with a 
standard sample. 

Anotner method, first used by Lanz sdorf 52} » 1s the pile 
oscillator. A mechanical drive oscillates a sample within the reactor. 
A neutron-sensitive ion chamber is used to record the overall power 


of the reactor. The oscillating sample causes an alternating 


eehementary definitions of the terms and concepts which are taken 
fron nuclear physics and reactor theory are given in Appendix I. 


“xumbers in brackets refer to the Bibliography. 








fluctuation in the overall reactor flux. The overall flux containing 
the variational component is sensed by the ion chamber. This signal 
current is treated so that the variational component is rectified and 
integrated with respect to time. The output voltae of the integrator 
is calibrated with a standard sample permitting the use of the device 
for measuring either cross sections or concentrations of impurities. 
In addition to the variation of overall flux caused by the 
oscillating absorber, there is a significant local variation. By 
placing the ion chamber close to the absorber the local signal can 
also be detected. If therral neutron cross sections are being 
measured, the best location for the oscillator and ion chamber is in 
a region such as the graphite reflector where neutrons are predominantly 
thermal. In Rue manner the overall reactivity changes induced by the 
sample are insignificant compared to the local disturbance. This 
technique was described by Hoover, et. al. {3} - The sample was 
oscillated through an annular ion chamber located in the graphite 
reflector of the X-1G reactor at Oak Ridge. The output of the ion 
chamber was amplified by a special feedback amplifier which passed 
orly the variational component of the overall signal. This output 
was fed to an integrator after anpropriate rectification. This method 
has been successfully used at Batelle Memorial Institute (4 } ° 
Comparing the pile oscillator method to the danger coefficient 
method, Langsdorf repurted an increase in sensitivity of ao factor 
of 2C over danger coefficient measurements with the same reactor. 
Using the Aerojet General Nucileonics 201 Reactor (AGN 201) for 
danger coefficient measurements, Carew and Mow { 5 made measurements 


as precisely as had been done with other reactors using pile 


oscillators. It therefore appeared possible that a pile oscillator 
for the AGN 2¢1 might produce even more precise measurements. 

It was decided to construct a pile oscillator for the AGN 201 
reactor and investigate its characteristics. Many possible designs 
had to be considered. In order to oscillate a sample through the 
core of the reactor, the glory hole would have to be used. The 
physical construction of the reactor would preclude placing a 
detector close enough to detect the local flux variation. If the 
oscillator were located in a region where a close-in detector could 
be placed, the access ports or thermal colusin would presumably be 
used. In this case the overall reactor flux would not be greatly 
affected and the reactor would be acting primarily as a flux field. 
All the previous devices of the type that oscillated within a thermal 
column or reflector were used in reactors having a thermal flux 
of the order of 10° neutrons pulenec Since the thermal flux of 
the AGN 201 is of the order of 10° neutrons ae nea this method 
would be at a disadvantage. 

Thus the type of oscillator which measured the changes in the 
overall reactor flux was indicated. From considerations of simplicity 
and reasonable expense, a design was base! on the employment of one 
of the permanent ion chninbers of the reactor and its associated 
amplifier on the Control console. It wns felt that the following 
procedure would be appropriate: 

1. To construct a mechanical drive mechanism which would 
oscillate a sample of material within the reactor with 
exactly reproducible motion over ench cycle and with 


accurite variable speed control. 





2e To devise an instruientation system which would reproduce 

accurately an input which sirulated the actual signal 

resulting from the oscillation of a typical sample. 

Se To perform experinents oscillating samples within the 

reactor. Using the instrumentation system devised, to perform 

measurements to determine if the output reading of the system 
is reproducible and to determine the sensitivity of the 


overall technique. 


2. Reactor kinetics for an oscillating absorver. 

The complete for:ulation of the effects of an oscilleting 
absorber ina reactor renvires extensive develon ert and involved 
expressions. i rigorous treatment 1s available frow Weinberg { 5} : 
However, since the pile oscillator method to be used with the AGN 20] 
reactor measures overall changes in neutron flux, the problem is 
consideranly simplified. Since interest is therefore confined to the 
ti: ewise variation of flux ratner t':an the spacewise variation, the 
case u.ny be considered in which the time-dependent part of the 
solution of the governing equations is separrxbvle fron the veo: etry— 
dependent parts. This is aciieved if the reactor is at or near the 
critical condition. 

In tie following treatment the assw ption of a homogeneous 
renctor is made. The one~—group space-independent kinetic enquations 


are: 


See Appendix I for develops.ent of these equations and definitions 
of syzbols. 





6 
dn Ak, {Oke an, _C. (1) 
dt Me g me roa 

os 


dc. wen 
aoe eo ee A; Cc. (2) 


In order to study the effect of an oscillating absorber, let it 

produce a change in k which is sinusoidal in form. Writing this as 
k(t) -Ak e!@* (5° = -1) 

and limiting Ok, to small values, A ki<« 1, the small amplitude 


variotions of n(t) and C . (t) are 


_ jJwt 
n(t) = n +n) e 
: n 
Cuties Ce + C: eae a. << 1 and n_ is the 
i 10 il no te) 


neutron density in the critica] state. Substituting these expressions 
and corresponding derivatives into equation (1) gives 
; 9 k , 
jeon, gwt _ Ak a Jot OA, n, e jot _aA on. _Jart_ e a 
® L © £ © £4 @ *%@ 
6 6 (3) 


-AX,P  ajet jert ) . jet 
ae - Jw JW 
n n, e + n,e + A; °io + A,o ie 


© @D i-l @ i=l ® 


If terms (3d and () are neglected as being secord order, and term (4) 


is neglected as deing small in comparison with term (2) » the equation 


reduces to 
6 


iW)n ~ “o Ak S 16} n. + ie (4) 
Saige j o 7 1 = 


Repeating the process for equation (2) gives 


qn 





jwt = x5 jit DS iE ae Jet 3 1 











JOC.) e = i] © + ~ ny + a no 
(5) 
a one jut 
ey eo + A, Oe ee 
Again neglecting second order terms gives 
ry, Bi 
Joti, = Ge >, 8, + iene A (6) 


Solving this equation for Ci and making the substitution into 


equation (4) gives 





6 6 
oe) ae > oo .\ ‘di Bi Alo Bo | (9) 
J pf oo gl Q| Jorn, Lj JwO+—”,; 
i= 1 ee 
Discarding the last term under the assumption Ak, n,<< n, and 
solving for ol . 
a) e 
O 
“ac 
no 8 (8) 
: eS cea Bay tn 
~ jo +d; 
i= 1 
Ak, 
- W(jw) Z 6 n 
| 1 Bi Ai (9) 
where W(jW) = i+ = | 3 = ag 
i=! 


Thus it has »Xeen shown, subject to the approximations made, 
that the ratio of the variational component of the neutron level 
to tne average neutron level is proportional to the input change in 
reactivity for a given value of G). W(jq) is called the reactor 


transfer function. 








The results of a set of calculations for different parameters 
4 , and for varying W is given in reference {7} - The particular 
results of interest and applicable to the AGN 2C1 reactor are shown 


in Figure lI. 
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4) (radians/sec) 


Amplitude of W(jw) as a Function of © 


Figure 1 





3. Description of pile oscillator for the AGN 2C1 reactor. 
(a) General arrangement. 
Figure 2 shows a block diagram of the pile oscillator 


constructed for the AGN 2C1 reactor. 


mac 


REACTOR 


CONTROL 
CONSOLE 





Figure 2 


Block Diagram of Pile Oscillator for AGN 201] 


The mechanical drive (1) moves the sample of material (2) 
through the core of the reactor with uniform linear motion. The ion 
chamber (3) detects the overall neutron flux level including the 
- chance induced hy the sawple and supplies a signal current to the 
linear amplifier (4) installed in the reactor control console. The 
output of the linear amplifier goes to the first stnges of the 
metering system (5) where the signal representing the response to the 
perturbation of the oscillating sarple is removed from the overall 
Signal. Tne positive part of this voltaze is integrated with 
respect to time by the intezrator (5) frou which the output may be 


read directly on a voltmeter or recorded on a strip chart recorder (7) 





for graphical analysis. The output voltage of the integrator is the 
criterion for the comparison of an unknown sample with a known sample, 
the frequency and number of cycles being fixed. 

(b) Mechanical components. 

In the design of the mechanical drive it was desired to provide 
an arrangenent which would be reasona}ly simple to manufacture and 
which would be adaptable to changes found necessary in the course of ~ 
the investigation. It was also desired to drive the sample in such a 
manner that the wave formr of the induced change in the multiplication 
factor would be sinusoidal. Since the absorption importance function 
of the reactor (see Appendix II) can be approximated by a cosine 
squared function, a linear wovement of the sample is equivalent to a 
Sinusoidal sarple motion through a linear flux gradient. Thus the 
mechanical drive was designed to produce linear movement of the 
sample through the reactor core. 

The drive mechanism is shown in Figure 3. A constant speed 
electric motor drives a hydraulic speed reducer to which is connected 
one sprocket for a roller chain. Another sprocket 1s mounted so as to 
stretch a heavy roller chain (3/4 in. pitch) between the sprockets. 
The connection of the chain is wade by a tool steel pin with a one 
inch overhang. This pin engages a vertical brass guide slot on an 
alurninum block which is free to move horizontally on two guide rods. 
Thus the travel of the pin causes the block to follow in a linear 
motion except at the end points where the direction of motion is 
reversed. The sample nolder is attached to the sliding block. 

The sample holder is made from 3/4 inch o.d. alurinw: tubing 


having a wall thickness of 1/32 inch. The tube was filled with 





Figure 3 


Mechanical Drive 
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reactor grade grapnite which was mac’iiined to close dianetral 
tolerances. A portion of the alu.inu. tube was cut away at the 
tidpoint. An interruption in the graphite at this point permitted 
the introduction of a test sample. The graphite filler was employed 
to .initize flux variations caused by the scattering of neutrons by 
the smuzple. 

(c) Instrumentation. 

The instrumentation systei converts the vnriations in neutron 
flux into useful infornation. The siznal is initiated in the ion 
chan ber in the reactor and amplified by the linear anrplifier on the 
reactor centrol console. Because of the low power level of the reactor, 
the ion chamber current at rated reactor power of one hundred 
milliwatts is ae of the order of nore anperes. The linear amplifier 
is used therefore to increase tne power level signal to usable levels. 
It is fron this point that the additional instrurentation takes its 
inform-tion. In speaking of the "signal" the reference will be made 
to the output of the linear amplifier. 

The signal voltage consists of the sum of three parts: a D.C. 
voltage representing the average reactor power level. an oscillating 
cowponent representing the response to the perturbation of the 
oscillating sample, nnd noise. Since the oscillating co. ponent is one 
or more orders of maynitude lower than the D.C. component, it. is 
necessary to seprirate it fron the overall signal to permit measurement. 
The frecuency of oscilixtion is qnite low naking removal of the 
variational component by usual amplifier techniques difficult. 

An anplifier of tine design used with the Oak Ridge oscillator {3} 


was constructed. This perforved well with low frenuencyv sinusoidal 


ll 





inputs but was not sufficiently stable when the input voltage 
contained a significant anount of noise. This was eventually 
abandoned. ; 

The probler was satisfactorily solved by employing operational 
amplifiers of the type employed in analog computer applications. A 
syste: was developed which subtracted the D.C. level fron the overall 
signal leaving the variational voltage. This voltage was rectified 
and integrated with respect to time. The integrator voltage after a 
selected number of oscillations was the basic value from the 
measurenent. 

In order to describe the instrumentation system in detail it 
will be necessary to first descride two elenentary analog computer 
circuits. 

A high gain operational amplifier is symbolized by , 


If the output is fed back through a resistor R, to the input and a 


t 


voltage e, is fed to the input through a resistor Ry: 


] 
R 
& 
Co 
the output voltage will be 

R R 
f : f 
aa R, el == G Qe where G = R, 


If adiitional voltages are fed into the input in like tanner, the 


negative sui Re R. will be formed. 
e rs comes + ae we * e@ * 
O Ry ] 9 2 


If R, = R, = Re etc, G is constant; the shortened symbol 
will be used for this case. 


If a capacitor of capacity C is used for feedbac: in place of 


l2 





the resistor, e e , the integral e — dt 


— 
~ 
>) 
oo 
) 
o-—— 
= 


will be formed. Likewise, additional inputs may be included 
R} | 
C 


oo AR - 


e@e¢ @ ] 
an oft e, + G, ey + | dt where G, = Rc 


If all the input resistances are evual, the shortened symbol > 


will be used. 

It was Buntea previously that the D.C. conponent of the signal 
veltage was to be found in order to subtract it from the overal] 
voltage leaving only the variational comriponent. This requires a 
definition of the variational voltage. This may be done in the 


e 
var 








following manners t 
| max 
Oo 
hat is, the tice integral of the variational voltage will not grow 


with time but will remain within the limit given. This may be 


expressed functionally: 








“ | 1 . I2, evar 
~e, — ~o a ie e,) aS RC max 
b RC), 


where e is the total signol, and en is the instantaneous value of the 
D.C. cowponent. In other words, if the integral of (e - e,) is formed, 
the value of the integral will not grow with time but will remain 


within the limit given. 


With this idea in mind the following scheme suggests itself 


oO iB) 
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Although the behavior of this system can be described in precise 
terminology, it ray be understood by considering the effect if i is 


not equal to the defined e If e | ls nunerically too small, there 


5° 
will be a net positive voltage input. This will generate a larger Py 
but, since this #8 negative by the sign reversing characteristic of 
the infezrator, the net positive input will be dirninished until 
equilibriun is established. 

hus far no account has been taken of the possibility that the 
average power level of the reactor may be drifting. This would have 
the effect of making e., drift. This can be minimized by selecting 
proper values of resistance and capacitance in the integrating 
anplifier so that 1t~ will follow the change in the average voltage. 
Generally speaking, the shorter the time constunt of the integrating 
aplifier, the better it will follow the drift. flowever, if the time 
const»nt is too short, it will tend to follow the oscillating 
component of the input as if it were a drifting D.C. Thus a 
compromise is reached based on the naxirium drift anticiprted,. 


With the D.C. component determined, the variational component 


can be found: 





In order to render the effect of drift insignificant, two more 
amplifiers was added repeating the above loop. The variational 
voltage is then rectified by a diode leaving only the positive halves 
which is the information signal, eC nfo’ The use of half wave 
rectification insures that the part of the signal which Nec ores 


e. corresponds to the sample position being in that part of the 


info 
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reactor where the cosine waveform is closely approximated (see 
Appendix II). The rectified voltage, Cfo? i8 fed as an input to an 


integrator. 


e. e. 
é info info 
y a r erent ee 3 oe 


| | | 


In order to check the reproducability of this system, a 


dt 


simulated input was fornied. This corsisted of an electronic function 
generator with a battery in series. This gave a D.C. voltage with a 
sinusoidal variation superinposed. The magnitudes and frequencies 
were chosen to correspond to the output of the linear amplifier. The 
system was allowed to reach equilibrium such that it balanced the D.C. 
corponent and then inteyrations were performed for a fixed number of 
cycles. This process was repeated, and oe resulting output voltages 
of the integrator were read. The readings were always the snme within 
a fraction of a percent. A slightly drifting voltage wes added with) 
no noticible change. 

A schenatic diagram of the amplifier system is given in Figure 4. 
Illustrations of the physical arrangement are shown in Figure 5. © 

A strip chart recorder was included in the instrumentation so 
that the finenl integrator output or any intermediary voltage could be 


é 


recorded. 
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Figure 5 


Instrumentation Equipment 
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4, Operation of the pile oscillator. 

In order to make a measurement, a standard specimen of iron was 
inserted in the sarple rod and the rod nositioned in the drive 
mechanism so that the sample was in the center of the reactor core 
when the drive mechanism was in its 1id-position. The amplifiers in 
the instrurent systev, were balanced and all inputs and condensers 
shorted. The drive mechanism was started ard the reactor control rods 
adjusted until the trace on the linear recorder on the control 
console showed no tendency to drift. This required several minutes. 
All shorts were removed from the instrumentation system except for the 
final integrating capacitor. After several minutes the instruments had 
steadied out and an integration could be cormnenced. This was done when 
the drive mechanism was in en extreme "out" position and continved for 
a fixed number of cycles after which the output of the integrator was 
recorded. The integration for the same sample was repented for 
approximately five readings. Then the sample was removed and replaced 
by a sample having an unknown percentage of impurity and the procedure 
repeated. In this manner the unknown sample could be relnted to the 
standard sample which was for all purposes free of the impurity. 

Several difficulties are inherent in the syster: which limit the 
total tine allowavle for the measurement of each sample. It is 
difficult to maintain the average level of the renctor at precise 
criticality. The nature of the instrumentation system is such that it 
must he balanced at least every half hour. The choice of component 
values in the instrumentation was to a larse part governed by the 
ability to compensate for drift effects. This in turn required that 
the oscillating frequency should not be too low or the instru ents 


would tend to follow tre oscillation as if it were D.C. drift rather 
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than the information bearing signal. 

The choice of oscillating frequency must take other factors 
into account. From the plot in Figure 1 it would appear that a lower 
frequency of oscillxtion is desirable. However, statistical 
considerations raise the possibility of advantage from higher 
frequencies since a greater number of observations can be accumulated 
in the same total tine. This arguzent assu.es that the relative error 
for each cycle does not increase as the frequency of oscillation is 
increased. Weighing all the contrasting factors it appears that a 
frecuency of oscillation near the inflection point of the curve in 
Figure 1 gives the best compromise. This corresponds to a period of 
about two seconds for each cycle. It should be noted that a complete 
cycle of the drive mechanism produces two cycles of reactivity 
oscillation. A set of tieasurements was made on the same sample at 
varying frequencies of oscillation but for the same totnl times. This 
confirmed the cycle time of about two seconds to be proper. 

In order to test the accuracy of the syster, a set of samples 
was chosen which had varying amounts of boron as an impurity in iron, 
These samples had been measured for impurity both by spectrographic 
methods and ay the danger coefficient method in the course of 
investigation reported »y Carew and Mow {5} - A typical group is 
tabulated in Table l. 


Table 1 
Samples of Iron Containing Boron 


Sample Weight (crams) Percentage of Boron’ Reactivity (+f ) 


A 6.381 0.0001 200 
B 5.351 0.6020 2106 
C 6.370 0.0050 232 
D 5.239 0.C072 218 
E 6.261 0.0250 267 
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9. Discussion of results. 

(a) Magnitude of error. 

In order to discuss the error of pile oscillator measurements, 
a sample standard deviation is employed as a basis for determining 
the sensitivity of the technique. The sample standard deviation will 
be here defined as the standard deviation of a single integration as 
deterviined frou several consecutive integrations with the same m 
material sample. The sample standard deviation can be projected to 
determine the sensitivity of a measurenent for which the integration 
process is repeated wit: the same material sample. The standard 
deviation of the arithmetic mean of the integrations would diminish 
by the square root of the number of repetitions. A set of typical 


readings for one sample is shown in Table 2 demonstrating the idea. 


Table 2 


Integrator Readings fron Integrations on Sample E 
Cycle time: 4.5C seconds; SO cycles each integration 


Reading (volts) 


29.3 
20.2 
28.C€ 
29.7 
25.9 


Average = 27.8 
Sample standard deviation = 1.7 or 6% 


Reactivity of sample 270“? approx. 
Maynitude of error 16¢P approx. 
” The error of a neasurerent would be diminished by 52 


approximately if five integrations were made on the same Sample. 
The result of numerous measurements such as the one tabulated was 
that the sample standard deviation for the method corresponded to 


an error of fron eee ey oc In order to reduce this »y a factor 


20 








of ten would require detem ining the arithmetic wen from about 
WC integrations for each sample. If each integration 1s as rapid as 
two minutes, over three hours would be required for the measurement 
on a single sample. External influences would undoubtedly come into 
play in’ this time. Employing the danger coefficient technique in a 
repetitive manner over a corresponding elapsed time with the AGN 201 
reactor gave a sensitivity of Ae 

Since tne instrumentation syster was capable of accurately 
reproducin., the test signal, exarnination of the experimental input 
sizsnal for variations is indicsted. The oscillating signal which is 
the input to the integrator was recorded on a strip chart recorder. 


An exact tracing fron a typical run is shown in Figure 6, 


| ji | 
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Figure 6 


Integrator Input Signal 


The severe variations are apparent. These variations are a 
consesjjuence of the total noise in the system. A discussion of this 
noise is desirable. 

(>) System noise. 

The total noise appearing at the final integrator input will be 
called the system noise. Noise is simply defined as unwanted signal { 8}. 
The system noise nay be traced to three principal sources: reactor 


noise, ion chamber noise, and amplifier noise. 


a | 








Altnourvh the equations of reactor kinetics are written for 
continuous functions, the actual physical processes occur in discrete 
steps. At sufficiently high power levels the nunber of events 
occuring are sufficient to limit deviations from the most probable 
pattern of events. When low powers or small perturdations are 
considered, random fluctuations becone significant. 

The reactor noise results from fluctuations in neutron flux 
which are caused by the discrete nature of the neutron creation 
process. It was noted that the delayed neutrons come from radioactive 
decay processes and are therefore subject to fluctuations. Both the 
concentration of delayed neutron precursors and their decay times are 
subject to randon variations. It is convenient to express the change 
in flux whic is caused by the oscillating sample in the forn of a 
comparison with the variations in flux not caused by the sample 
(i.e. reactor noise). This is simply a signal to noise ratio. 

Frisch and Littler {0} derived an expression for the ratio of 
signal to noise for a single observation of an oscillating sample: 


is 


27N, T A (9) 4 





S(@) = (Alc), 9 where 
1+TA (w) 


(Ak) - = change in multiplication produced by the sample 
vT = probability of the detector detecting a fission 


t 


mean fission rate 


N 


f 
T = total tire of the observation 
ac. ) = amplitude response of the reactor. 
Rewriting this as S(w) = (Ak) 2N, T G —s where 
ae A(@) i 
Tos te) 


shows the signal to noise ratio with respect to renctor noise 


ae 





improves as the square root of the ratio between two power levels. 

Ion chanber noise has been demonstrated experimentally {10} to 
vary inversely as the square root of the ion chamber current. Since 
current 1s proportional to neutron flux, this varies in the same 
manner as the reactor noise. Thus the total noise appearing at the 
linear amplifier varies inversely as the square root of the reactor 
power level. 

The third type of noise is that produced in the electronic 
equipment. This was found to be printarily sixty cycle hum. This was 
attenuated in the measuring scheme employed. 

The precision of a measurement using this technique will be 
es-entially determined by the system noise described. Since this 
increases as the power level of a reactor is reduced, it will be a 
formidable obstacle to the success of a pile oscillator on a low 
power reactor such as the AGN 2Cl. 

In reference {11} an expression was derived which is similar to 
that of Frisch and Littler discussed previously. In this case, 
account was taken of the correlation involved in the type of 


experiment. The expression for signa] to noise ratio is 





S = a 2F with F the mean fission rate of the 


N x, 


reactor, T the total time of the observntion on the given sample, 
and Xa factor involving the mean number of neutrons released per 
fission. For a typical nieasurewent (a single integration), this 


expression gives: 


(.l watts x 3 x io! fissions/watt sec. 


F = 
= 3x 10° fissions/sec. 
T=2.0 x 107 sec. (a typical value) 
eS 
= 2 10° Ak 
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Thus for a signal to noise ratio of unity, k must be of the 
order of ome which 1s approximately ten micro-rho. It would anpear 
that the linitation on the results has been placed largely by the 
reactor noise. 

The effect of the noise on the danger coefficient method is slight. 
In the course of an observation, the noise is involved only when it is 
necessary to deternine the reactor period. If the criterion for this 
determination is the position at which the trace of a linear amplifier 
output on a strip chart recorder crosses a given ordinnte, the average 
value of the flux level can ve rend through the noise by drawing in a 
line which is the best estimate of the average signal and accurately 
o»serving where this line crosses the ordinate. A human operator has 


considerable ability in finding the average value of such a signal 


through noise. This technique is illustrated in Figure 7. 
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Intersection accurately determined 
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Figure 7 


Time Determination from Strip Chart Recorder Trace 


Thus the statistical advantage of ao pile oscillator over danger 
coefficient observations is open to considerable doubt. The noise 


factor may be so disturding as to outweigh all other considerations. 





6. Conclusions. 

It is concluded that a pile oscillator for the AGN 2C1 reactor 
utilizing the installed ion chamber and linear amplifier does not 
give satisfactory results. It is not possible to attain with the 
equipment employed the precision oytainable with danger coefficient 
measurements on the same reactor (and using less eouiprent). 

Since the instrumentation system consistently reproduced a test 
input, output variations with the experiwental input are attributed 
to the nature of the information contained in the input siznal. The 
reactor noise appears to de the major cause of the variations. 

In order to improve the sensitivity of the pile oscillator, 
either one or both of two things must be done. The quality of the 
input siynal could be improved vy the use of a better ion chamber 
placed closer to the core. Placerxent close to the path of oscillation 
might pernit detection of the local signal to incrense signal to 
noise ratio, Another inprovement would result from operating the 


reactor at a higher power level than the presently licensed limit. 
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APPENDIX I 
ELEMENTARY BACKGROUND 

In a nuclear reactor, the reaction senuence »y which neutrons 
are produced and lost is self sustaining if proper conditions exist. 
This is described by the multiplication factor, k, being unity. This 
is the "critical" condition for the reactor. If k is greater than 
unity the reactor is supercritical and the neutron population will 
“increase. Conversely, if ranils less than unity, the reactor is 
subcritical and the neutron population will decrease. 

Physical materials have properties which are important in 
relation to reactor phenomena. The probability that a material 
will interact with specific nuclear particles is descrivded 
quantitatively in terms of the "cross section" for the reaction. 
Of particular interest is the probability that an atom of the 
material will absorb a neutron. This is expressed in terms of 
the microscopic absorption cross section and denoted by ga . It 


2 2 
= cm. The product of 


carries the units of "barns"; one "barn" is 10 
the number of atoms in a unit volume and the microscopic cross 
section is terned the tiacroscopic cross section and denoted by a ° 
Thus eS = NG, - N is usually expressed in units of atoms per cubic 
centimeter; therefore 2. has the units of we The number of 
neutrons in a given volume is termed the neutron density and denoted 
by n. It is sufficiently accurate to state that 1f the total neutron 
population is in equilibrium, the neutron density (neutrons per cubic 
centimeter) is similarly in ecuilibrium. Thus on = 0. The partial 


derivative is used since n is also a function of space coordinates. 


However, for a reactor at or near critical, the ceoretrical shape of 





the distribution of n will not change significantly. 

Since the energy or velocity of a neutron is important in 
describing its behavior, the product nv where v is the neutron 
velocity is useful. The term neutron flux is defined as this product 
and is of necessity a function of the energy (or velocity) considered. 
This is denoted by the symbol § and carries the units of neutrons 
oe a If the neutrons are in thernal equilibrium with the 
surrounding medium, they are called "thermal neutrons" and the flux 
is called "thermal flux". Since thermal neutrons cause most of the 
fissions, knowledge of cross sections at thermal energies is of 
considerable importance. Cross section data are usually given for 
neutrons having a velocity of 22C0 meters per second, the most 
probable velocity per unit velocity for thernal neutrons at room 
temperature. 

If a reactor is not critical, there will be a net change in the 
nucber of neutrons. This may be traced to a corbinnation of three 
causes: leakage out of the reactor, absorption within the reactor, and 
a gain fro. sources of neutrons. The leakage per unit volu.e is 
ciaracterized by a term -DV ¢ arising from diffusion theory; D is the 
diffusion constant. The absorption term is written 2 if utilizing 
the definition of macroscopic absorption cross section. The source 


term may be written sinply as S. The resulting equation for the 


neutron balance per unit volume per second decones 


co pV -ZS f+s (1) 


This is known as the time-dependent neutron diffusion equation. 
All but a small fraction of neutrons fro. fission are given off 


immediately with the fissioning of urnaniur nuclei. These are terned 








prompt neutrons. The small remaining fraction (C.0065 of the total) 
are not emitted inmediately but follow the radioactive decry of 
certain fission product nuclei. A total of six species of fission 
product nuclei which are precursors for the emission of these delayed 
neutrons are known for ues Of the total neutrons emitted per unit 
tire, the fraction which is due to eaci species of delayed neutron 
6 

precursor is denoted by 6 i? and GC = Da fe — 0.0065. 

The remaining possible sources of neutrons are from external 
means such as a radium-beryllium source or cosmic rays. 

The total source term way be written: S = =F + Sa t Sy, for 
prompt, delayed, and primary or external neutrons respectively. 


The important characteristics of the delayed neutron species or 


groups are given in Tavle 3. 


Table 3 


Delayed Neutron Characteristics 1 


Group Index (i) Decay Constant, A, (wee *) Fractional Yield, B, 





] 0.0127 C.C00247 
2 0.6317 0.00138 
3 C.115 ©.€0122 
4 0.31] 0 .€0265 
5 1.40 C.0( C832 
6 387 C.C00169 


Total: @= 0.00650 ( =236,) 


It is important to note that the decay is the inverse of the 
averave lifetime and that the actual emission of tle neutron can 
occur over a time distributed about the average value. Also the 


fractional yield of a given species is an average and o]s0 is 
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suodject to fluctuations. 

The discussion of the time behavior of a reactor is grently 
simplified if the materiirls in the reactor core are homogeneous. This 
is approximately the crse with the AGN 201; therefore this assumption 
is used. 

The technique of the pile oscillntor depends on the reasurerent 
of the change of the reactor rultiplication f»xetor resulting from the 
oscillating oabsorver of neutrons. The effective multiplication factor 
of the reactor, k, is defined as the ratio of the number of neutrons 
present in the renctor at a given time to the nurber of neutrons 
present in the immediately previous generition. The concept of generation 
implies a characteristic lifetime for the neutrons which is the average 
time between the "dirth" and "death" of a neutron. Since neutrons can 
be lost fror. the reactor either ay absorption in sone nuclei or by 
leakage frou the reactor boundaries, it is useful to separate the 
effects. Thus a lifetine, f,, can be postulated for thermal neutrons 
in an infinite syste:. from which there can be no leakage. The 
reciprocal of the macroscopic absorption cross section is equal to the 
mean free path of the neutron; that is, the average distance which the 
neutron travels until it 1s avsorbeil. Since >_ has the diuensions of 

a 
ew, the mean free path has the proper dimensions of centimeters. If 


v is the characteristic velocity of a thermal neutron, the neutron 


lifetine Q,is then —- - If a factor P 


Ve t 


probability that a thermal neutron will not leak from the reactor is 


which represents the 


introduced, the lifetiie for finite geonetry is written f = p Pie 
=) 
This concept assures that all the neutrons diffusing in the reactor 


are thermal neutrons. A more correct expression would include the 
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time a neutron spends slowing down from fission energy to thermal 
energy. Since this time is one or more orders of maynitude smaller 
than the thermal lifetime, its neglect will not he overly serious. 


The factor k is related to the multiplication factor k for 


Qo 
an infinite reactor of the sa:re core material vy employing an 


additional factor P, correcting for the leakage of neutrons durin 


f 


the slowing down process. Thus k = k,P,P 


o PP, The expressions relating 


these factors to the materinl and geonetry of the reactor are not 
nertinent to this treatment. Fron the definition of k it is seen that 

if k = 1 the reactor is critical; that is, the rate of creation of 
neutrons is equal to the rate of loss of neutrons and the time variation 
of neutron density is zero. The cuantity Ak =k - 1 is a useful 
measure of the tire behavior of the neutron density. If Ak is not 

zero, the reactor will not be in a steady state condition. The time 


behavior of the reactor may also be related to the term "reactivity" 


At 


7a A convenient unit is the micro-—rho 


defined by the relation @ = 
where ae = ¢ Tor >. If k is very close to unity, ie = eC. 

In a finite reactor, if the neutron density at a given time is n, 
there will be Im neutrons after the lifetime / - Thus there will he a 


net increase of Akn, or a production rate of = n if all the 


neutrons are pro:pt. To correct for the delayed neutrons the factor 
k A 


must be sudtracted. Finally there will be added those neutrons which 


ny ere ( is the total fraction of neutrons which are delayed, 


come from the decay of the delayed neutron precursors. If C. is the 
concentration of a given delonyed precursor and Se is its 

© 
characteristic decay corstant, the rddition will be A, C. 


. a= 
where the sun. is over the six known delayed neutron groups. 
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The eruntion for the neutron bnlance is thus 


dn Sk k 6 v C (2 
cs 7 = ee i >) 


A=! 





This equation involves the unknowns n and C.. A second ecuation may 


be written for the fornation of the delayed neutron precursors: 


ac 
1 k ; 
rae os e 1 n a Ts C; (3) 


J 


Equations (2) and (3) are the equations of the so called one-croup 
space-independent rmiodel of reactor kinetics. 

A development starting from tne diffusion equation and using a 
more sophisticated analysis {12} which takes some account of the 
slowing down of neutrons from fission energies to thermal energies 


wives the ecuation: 
(8) y | =) elu 
(i(1 B) 1} nv a A, c, = f * (4) 


The factor p is the fraction of neutrons which escape absorption 
while slowing down from fission to thermal enervies. Using the 


definitions: keky Pi P. 





f = f,P4 equation (4) becomes 
aa -t f n ~ L n + 4 i 


/ 
where C. = p Pe Ci. 


Thus ecuation (2) will be in asreenent if the C, are defined as the 


teffective" concentrations of delayed neutron precursors. 
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APPENDIX II 
DESCRIPTION OF THE AGN 201 REACTOR 
, The AGN 201 is a scall low power reactor having a licensed 
output of one hundred milliwatts. With incrensed shielding this can 
be increxsed to twenty watts. The reneral chnracteristics of the 


reactor are given in Table 4. 


Table 4 


Characteristics of the AGN 2C1 


Power (licensed) 1CC6 milliwatts 
haximunci Thermal Flux @€ 10C mw 4.3 x 10° we sec a, 
Critical Nass 665 em Cae (20% 
enrichment in UO0,) 

Core Size (cylinder) 10" dia. x 94" nen 
Moderator (polyethylene) 11 ke 
Reflector (graphite) 20 em thick 
Shield (lead) 10 em thick 

(water) 55 em thick (1000 gals) 
Overall Diameter 6.5 ft. 
Overall ileigcht SI i ar 


The “clory hole" is a one-inch aluminum tube which runs 
horizontally through the entire reactor and passes through the 
center of the core. This permits introduction of various materials 
into the region of highest neutron flux. A cross section view of 
the rexctor along a line through the glory hole is shown in Figure 8. 

The variation of thermal neutron flux with radinl pcesition in 


the glory hole is shown in Figure 9, 


33 








Graphite 


Glory Hole 


Figure 8 


Reactor Cross Section 


View 
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Figure 9 


Therval Nevtron Flux Through Glory Hole 
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The rate of absorption of themial neutrons by a material! is 
proportional to the ther al neutron flux throuch the material. The 
therral neutron flux varies with position through the reactor. Since 
k is proportional to the rate of absorption of neutrons, k will 
depend on the position of an absorber in the reactor. The relative 
effect of an absorber can be referred to its effect at the center of 
the reactor. This is expressed as an importance function. The 
variation of absorption importance function with position in the 
glory hole is shown in Figure 10. A curve of squared cosine function 
is fitted to the importance function curve over the dianeter of the 


core. 
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Absorption Importance Function 
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Figure 10 
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